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PACS 68 . 15 . +e - Liquid thin films 

PACS 68 . 03 . Cd - Surface tension and related phenomena 
PACS 68 . 03 . -g - Gas-liquid and vacuum-liquid interfaces 

Abstract. - Within the framework of a semi- microscopic interface displacement model we analyze 
the linear stability of sessile ridges and drops of a non-volatile liquid on a homogeneous, partially 
wet substrate, for both signs and arbitrary amplitudes of the three-phase contact line tension. 
Focusing on perturbations which correspond to deformations of the three-phase contact line, we 
find that drops are generally stable while ridges are subject only to the long-wavelength Rayleigh- 
Plateau instability leading to a breakup into droplets, in contrast to the predictions of capillary 
models which take line tension into account. We argue that the short-wavelength instabilities 
predicted within the framework of the latter macroscopic capillary theory occur outside its range 
of validity and thus are spurious. 



Introduction. — Many practically relevant processes 
involve liquids covering solid substrates. The liquid may 
be a paint, a lubricant, an ink, or a dye. The solid may 
be hard as a metal or soft as human skin or gel; it may 
also show a homogeneous surface or be finely divided, as 
it is the case for suspensions, porous media, or fibers. Due 
to the complexity of these systems even extensive research 
in the past has still left many open questions concerning 
in particular the intricate phenomena taking place in the 
proximity of contact lines formed at three-phase bound- 
aries [1-4]. 

Considering a small liquid droplet deposited on a flat 
solid surface, two distinct types of behavior can be ob- 
served depending on the sign of the so-called spreading 
coefficient S ~ o~sg — a sl — c, where <jsg, &SL-, and a 
are the tensions of the three interphases which meet at 
the contact line: solid-gas, solid-liquid, and liquid-gas, re- 
spectively. If, for a nonvolatile liquid, S > 0, the liquid 
tends to shield the solid from the vapor phase which cor- 
responds to the complete wetting situation; the drop will 
thus spread, tending to cover as much of the solid surface 



as possible. Conversely, S < corresponds to the par- 
tial wetting regime; the drop will then evolve at some rate 
(see, e.g., ref. [5] and refs. therein) towards an equilib- 
rium shape characterized by a local thickness profile h(r) 
as function of the lateral coordinates r and an equilibrium 
contact angle 6>o at which the liquid-gas interface meets 
the solid surface. If the drop is small enough for gravity 
effects to be negligible, it assumes the shape of an ideal 
spherical cap. Typical examples of partial wetting situ- 
ations are droplets of dew on leaves or water drops on 
greasy glass. 

Two centuries ago Young [6] considered the balance of 
local surface forces at the edge of such a drop and deduced 
that the equilibrium contact angle 9y must satisfy 



costly = 



&SG — &SL 

a 



(1) 



(a) 
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For a long time, eq. (1) served as a tool for inferring surface 
tensions from measurements of contact angles By . On the 
other hand, it was also understood long ago that Young's 
equation is strictly valid only for the somewhat idealized 
geometry of a liquid "wedge" - a situation in which bulk 
phases are of macroscopic extent, the interfaces are flat, 
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and the contact line is straight. Thus, for axisymmetric 
drops of base radius R, a modified Young's equation holds, 
which takes into account, in leading order, the finite drop 
size: 

T 

cos 9 = cos 9y 5 , (2) 



where r is the excess free energy of the solid-liquid- vapor 
system per unit length of the contact line [7,8]. As the one- 
dimensional analogue of the surface tension, r is usually 
referred to as the contact line (CL) tension. 

Formally, 9 in eq. (2) reduces to By given by eq. (1) 
only in the limit R = 00 (straight CL) or if r = 0. Typical 
values of r are of the order of 10~ 10 to 10~ n J/m (see, e.g., 
refs. [9-14] and refs. therein) and may entail only small 
deviations of 9$ from 9y for generic values of R. However, 
there might be other interesting phenomena caused by the 
very fact that r is nonzero. We note that whereas surface 
tensions are necessarily positive, r is not constrained to 
be positive and may adopt negative values as well [12,15]. 

The most obvious consequence of a negative CL tension 
is that sufficiently small drops no longer tend to coalesce, 
as soon as the negative CL tension overcompensates the 
ensuing gain of surface free energy. This critical drop size 
is of the order of r/er, which is typically in the nanometric 
range. Over the past few years, this nanoscale has become 
directly accessible to experimental observation, which has 
renewed the interest in negative CL tensions and their 
implications [8,10,16-22]. 

A central issue of these studies is that of CL stability. 
An early analysis has been performed by Dobbs [16] within 
the framework of the interface displacement model (see, 
e.g., ref. [17]). Considering the ideal wedge geometry, he 
has shown that a negative CL tension does not entail a CL 
instability with respect to small amplitude deformations 
of any wavelength. 

The question of the CL stability has been further ad- 
dressed for other, experimentally more relevant geometries 
such as axisymmetric drops [18] or long liquid filaments 
(ridges) [19,20], within the framework of a macroscopic 
approach based on the second variation of the free en- 
ergy composed of the interface and line contributions [21]. 
The latter authors have obtained complete phase diagrams 
as function of #0, t, and the perturbation wavelength A, 
showing both stable and unstable regimes. As opposed 
to the interface displacement model, the macroscopic ap- 
proach with negative CL tension yields a range of unstable 
short wavelength modes, bound from above by the critical 
wavelength 

a sin 9o 

Remarkably in ref. [20] a negative r was found to desta- 
bilize a liquid filament with respect to contact line defor- 
mations of any wavelength for a sufficiently small contact 
angle 0q. As it stands this result queries the standard view 
that the spontaneous break-up of a filament into a chain 
of droplets (known as varicose or Ray leigh- Plateau insta- 
bility) occurs only for wavelengths large compared with 



the lateral size of the filament. Additionally, in ref. [20] 
a positive CL tension was shown to stabilize the varicose 
in accordance with the observation that the critical (i.e., 
marginally stable) wavelength increases with r. 

We note, however, that the analyses in refs. [18-20] treat 
a, 9q, and r as independent parameters, whereas for an ac- 
tual system they are linked on the microscopic level. Con- 
sequently, realistic partial wetting situations correspond to 
a certain subspace of the surprisingly rich phase diagrams 
obtained by the authors, and one may expect a micro- 
scopic or semi-microscopic analysis to clarify the picture 
provided by those macroscopic models. 

This so far unnoticed contradictory state of the litera- 
ture prompted us to revisit the problem of CL stability for 
liquid ridges and drops within the framework of the inter- 
face displacement model. Our aim is to determine what 
kind of instabilities can be encountered in realistic models 
and hence, which parts of the phase diagrams deduced in 
refs. [18-20] are physically plausible. 

In the following we first introduce the interface displace- 
ment model, together with general considerations concern- 
ing the stability analysis. We then focus on specific geome- 
tries, i.e., ridges and drops, analyzing both with respect 
to perturbations leading to a deformation of the contact 
line. We summarize our results in a final section. 

Basic properties of the interface displacement 
model. — We consider a two-dimensional, homogeneous, 
and flat substrate completely covered by a thin liquid film; 
h = h(x,y) describes the local film thickness (film in- 
terface displacement relative to the substrate) at a point 
(x, y). $(/i) denotes the effective interface potential, which 
incorporates the net effect of the substrate-liquid and 
liquid- liquid interactions [1,23]. We consider a general 
form of such potentials which might exhibit several min- 
ima and maxima; we only stipulate that $(/i) vanishes for 
h — > 00, is differentiable, is bound from below, and has its 
deepest minimum at a microscopic thickness h = a > so 
that $'(a) = 0. A typical example of such a potential is 



$(/>) = 



A 



12nh 2 



1 



1 ra 
4 \h 



(4) 



which corresponds to Lennard- Jones pair interactions, 
disregarding other correction terms to the leading term 
<~ h~ 2 [23]. A < implies partial wetting. 

The free energy per surface area of a planar liquid film 
with prescribed thickness h is given by 0(/i) = <jsl + cr + 
3>(/i) so that (7 sg — minfl(h) = ctsl + c + 3>(<x), i.e., 

h 

<f> (a) — S. Accordingly, $(a) < corresponds to partial 
wetting. 

This description formally applies only to films of con- 
stant thickness. For laterally varying thickness profiles 
h (x, y) the model relies on the assumption that the two- 
dimensional gradient V = xd x + yd y (i.e., the slope) of 
the profile h is small, i.e., |V/i| <C 1. In this case the flat 
film expression of the interface potential &(h) is locally 
valid. 
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The effective interface Hamiltonian £ of such a smoothly 
varying thin film has the form [17] 



£ [h] 



dxdy 



\a(Vhf 



<f>(h) 



(5) 



in which the surface tension a reduces to a scaling parame- 
ter for A is the surface area. An important assumption 
is that the liquid reaches the whole substrate, so that the 
nonwet surface is covered by a film with a microscopic but 
finite thickness a. 

The stationary profile h(x,y) minimizes the functional 
in eq. (5) under geometry-specific constraints. As in any 
meso- or microscopic approach, there is actually a con- 
tact line region within which the interface interpolates 
smoothly between the characteristic features of h(x,y). 
In the case of the liquid wedge geometry (invariant along 
the y axis), the asymptotes of h(x) intersect such that the 
contact angle 6q (= 9y in this case) is formed [22]: 



6> = arccos [1 + $ (a) /a] 
9 a w v/-2$(a) /a. 



(6a) 
(6b) 



The expressions in eq. (6), with the latter being valid 
for 6>o <C 1, are consistent with the general expression of 
Young's law. Furthermore, the CL tension r of a wedge is 
given by the de Feijter and Vrij formula, referring to the 
excess tangential force acting on the contact line region [9] 
or the excess energy localised in that region [22]: 



OO 



dh 

dx 



dh 



= V2^ / V*(/i)-*(a) - y/-$(a) dh . (7) 



Here h{x) minimizes £ [eq. (5)] with the boundary condi- 
tions d x h (x — > — oo) = and d x h (x — ► +oo) = 8 a . 

As pointed out in refs. [9,22], the interface displace- 
ment model allows for arbitrary values of the line tension 
[eq. (7)] by means of appropriately chosen <E> (h) (slow de- 
cay, large humps, etc). However, in particular for nega- 
tive line tensions, there is a close relationship between the 
contact angle 9q, the line tension r, and a characteristic 

thickness £ = ^ / a + °° [#o - aH x=R -i (fc) ] d ^ , related to 
the range over which the interface potential $ influences 
the shape of the contact line. As an example, for the effec- 
tive interface potential in eq. (4) r is negative and £q ~ a - 
In terms of this effective range £q one has r = —<t9q£q 
and the critical wavelength Ao of CL instabilities as pre- 
dicted by eq. (3) (see, e.g., ref. [20]) is given by Ao = £o/0o, 
i.e., it is comparable with the lateral width of the contact 
line region within which the shape of the liquid- vapor in- 
terface is influenced by <&. Thus the capillary model eq. (3) 
is based on predicts an instability at a length scale for 
which its use is rather questionable. This motivates and 



encourages us to investigate the stability of ridges and 
droplets by using a more microscopic model which is able 
to cover the range of wavelengths below Iq/Oq. 

Stability analysis within the interface displace- 
ment model. — For the effective Hamiltonian in eq. (5) 
the energy variation £ = £ [h + e^~\ — £ \h\ due to a small 
perturbation e^{x,y) around a stationary profile h(x,y) 
is given, to second order in e -C 1, by the expression 

£ = e jj dxdy [$' (h) + a (Vh) • (V*)] 



dxdy 



$" (h) * 2 + a (V*) 5 



(8) 



with $'(/i) = d$/d/i and = d 2 $/d/i 2 . The require- 

ment that the first variation vanishes for any perturbation 
'J preserving the volume of the liquid, i.e., with the con- 
straint 

(9) 



(10) 



■^/dxdy = , 
yields the Euler-Lagrange equation 

$' (h) - <rV 2 h = P 

for the equilibrium profile h. The Lagrange parameter P 
fixes the volume and plays the role of the Laplace pressure 
so that for the homogeneous equilibrium film thickness 
h(r) — a one has P = 0. 

The stability condition is that the second variation 

£ = ie 2 J J dxdy [$" (h) * 2 + a (V*) 2 ] (11) 

is non-negative for any perturbation '5 fulfilling eq. (9). 

For a suitably normalized $ , the integral on the right- 
hand side of eq. (11) is the expectation value of the Her- 
mitian operator 



U(x,y)-aV 2 



(12) 



where U (x, y) = <&" [h(x, y) \ . We call P the "pressure op- 
erator" because eP^> is the additional local pressure gen- 
erated by the perturbation e'S (c/ eq. (10) for stationary 
contributions). 

Analogous to quantum mechanics, \& can be regarded 
as the "wave function" of a particle and P as its "Hamil- 
tonian" with potential energy U (x, y) and kinetic energy 
(— crV 2 ). Thus the linear stability analysis of the station- 
ary profile h reduces to an eigenvalue problem for P: 



(13) 



which has the structure of the time-independent 
Schrodinger equation. As it will turn out, the elemen- 
tary perturbations of interest automatically preserve the 
volume of ridges and drops, and thus in every specific case 
the constraint in eq. (9) is satisfied. 
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The eigenfunctions with positive energy E are called 
stable, while the ones with E < contribute to the in- 
stability of the system. Thus the stationary solution h is 
only stable if P has no eigenfunctions with E < 0. 

A basic feature of the Schrodinger equation is that its 
ground state has no zeros inside its domain of definition 
(see § 20 in ref. [24] and chap. 6 in ref. [25]). This holds 
not only for the 2D Hamiltonian P, but also for the effec- 
tive ID Hamiltonian(s) it reduces to for geometries with 
translational or rotational symmetry. 

Another basic feature, valid for any stationary film 
structure on a homogeneous substrate, is the indifference 
with respect to lateral shifts in any direction u of the (x, y) 
plane, following from the translational symmetry on a ho- 
mogeneous substrate ($ does not depend on x or y explic- 
itly). Specifically, for any unit vector u, * = u • Vh is an 
eigenfunction of P with eigenvalue 0: 

P (u • Vh) = $" (h) (u • Vh) - aV 2 (u • Vh) 

= u • V [*' (h) - aV 2 h] = 0. (14) 

S v ' 

P— const 

The knowledge of the explicit expression W = u • Vh for 
lateral shift modes with zero energy turns out to be very 
useful whenever the stationary solutions exhibit specific 
geometrical properties (e.g., ridges or drops). This forms 
the core of our following analysis, together with exploiting 
the existence of a unique ground state with no zeros. 

Ridge stability: preliminaries. — If the liquid 
structure is translationally invariant along the y axis, the 
stationary solutions are of the form h(x) and eq. (10) re- 
duces to 

&(h)-a8%h = P. (15) 

The following analysis is based on solutions h(x) of eq. (15) 
which resemble a liquid ridge centered symmetrically 
around x = with h(\x\ — > oo) = ho and d x h (|x| — > oo) = 
so that P = ®'(ho). Thus the first integral of eq. (15) is 

\a (djif (h Q ) -(h-h )P. (16) 

A typical solution of eq. (16) is shown in fig. 1 (solid curve). 

In view of studying eqs. (12) and (13) in the present 
case for which U(x,y) = U(x), we can restrict the linear 
stability analysis to the Fourier components of the defor- 
mations: 

*(x, y) = V>(x; fc) cos [fc y + a(k)] , k ^ , (17) 

with arbitrary phase shifts a(k). For k ^ 0, W automati- 
cally satisfies the volume conservation constraint in eq. (9) . 
Thus for k ^ eq. (13) reduces to the ID Schrodinger 
equation: 

Hip (x) = Eip(x) , H = H + ok 2 , (18) 
H = U(x)-ad 2 x , U (x) = $" [h (x)] . (19) 



If -0 is an eigenfunction of H with eigenvalue E, it is also an 
eigenfunction of H with eigenvalue E = E + ak 2 . Due to 
a > this relation allows us to infer the stability from the 
limiting case k — > 0, which is the least stable mode. Thus 
the 2D mode \I/(x, y) = ip(x) cos (fc y) is stable (E > 0) for 
all fc, unless E < 0, in which case there is an instability 
for |fc| < (-E/af 2 . 

With that at this point the general results are: 

• Within the framework of the interface displacement 
model a stationary solution of eq. (15) can be unstable 
only with respect to long-wavelength perturbations 
(such as Rayleigh-Plateau instabilities). 

• The discussion of unstable modes can be reduced to 
that of eigenfunctions ip of H with eigenvalues E < 0. 

This reasoning, which does not depend on the form of 
Q(h), rules out linear instabilities of ridges with respect 
to short-wavelength deformations of the contact line. 

The following section provides a more detailed study 
of the long-wavelength instability of ridges, including the 
limiting case of a wedge, which can be viewed as a ridge 
of macroscopic transverse size. Finally we shall turn our 
attention to axisymmetric drops. 

Ridge stability: detailed analysis. Figure 1 
shows the vertical cross-section along x of a ridge (solid 
curve) centered at x = and aligned with the y axis. The 
basic features are an apex at x = and a symmetric de- 
cay on both sides towards h(— oo) = h(+oo) = ho. The 
thickness ho corresponds to a stable film, i.e., $" (ho) > 
so that U (±oo ) > 0. For eq. (4) this implies ho < v^3a. 

By analogy with quantum mechanics, an eigenstate ip 
of H with eigenvalue E < is bound by the potential 
U(x) = $" [h (x)] . According to the so-called "oscillation 
theorem" (see, e.g., § 21 in ref. [24]), the bound states 
are non-degenerate and indexed by their number of nodes. 
Since the stationary solution h(x) and thus the potential 
energy U(x) are symmetric with respect to x = 0, the 
eigenstates are either even (symmetric) or odd (antisym- 
metric). 

As a special case of eq. (14), ipi (x) = d x h(x) is an 
eigenfunction of H with eigenvalue E\ = 0. The asso- 
ciated 2D deformation [eq. (17)] corresponds, for fc = 0, 
to a lateral shift of the ridge, and for fc ^ 0, to finite- 
wavelength meandering. Since is an eigenmode of H 
with eigenvalue E\ = ak 2 , a ridge is generally stable with 
respect to meandering. 

The eigenfunction if>i has a single node at x = 0, i.e., 
at the apex of the ridge. From the properties of bound 
states we infer a unique ground state ipo for H with an 
energy Eq < E\. The mode ipo has no nodes: the associ- 
ated 2D perturbation corresponds, for fc ^ 0, to periodic 
bulging and shrinking of the cross-section of the ridge, and 

1/2 

is effectively unstable for |fc| < (— Eo/a) (varicose or 
Rayleigh-Plateau instability). Perturbations with wave- 

— 1/2 

lengths shorter than A = 2ir/ (—Eo/a) are stable. As 
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Figure 1: Stationary profile h (x) /h (0) of a liquid ridge (with 
h (0) ~ 2.81a and ho = 1.04 a) and corresponding potential 
energy U (x) = $" [h(x)] (eq. (12)) of the ID Hamiltonian, 
normalized by U(oo) ~ 4.92a9 2 /a 2 . Also shown are the two 
eigenstates of P of lowest energy, ipo and tpi, normalised by 
their maximal values, h, U, and ipo are symmetric and ipi is 
antisymmetric. 

the only eigenmode of H with E < 0, for any ridge Vo 
is the only source of instability. Figure 1 shows both ip\ 
and tpo for a ridge, together with the stationary profile h 
and the symmetric potential well U (x) — <I>" \h (a;)] . <& is 
given by eq. (4), i.e., the line tension is negative. 

Unlike for ipi , there are no general analytic expressions 
for ipo and its energy Eq. Reference [26] considers the 
limit of large ridge width L x — |4a/6* and small internal 
pressure P = 5$ (a) /a, where 5 is a small dimensionless 
parameter. In that limit, it can be shown [27] that the 
marginally stable wavelength Ao scales linearly with the 
ridge width L x . The key result is that the limiting ratio 
lim Xq/L ~ 2.62 is intrinsic: it depends neither on details 

L — >oo 

of $ (h) nor on the value of r as determined by eq. (7). 
Interestingly, ref. [20] also analyzes the limit 9o — > and 
reports A /L — 2.62 for the specific case in which the 
CL tension, used as an input variable in that macroscopic 
approach, is exactly zero. 

The linear scaling of the marginally stable wavelength 
with the ridge width is consistent with the marginal sta- 
bility of the wedge, viewed as a ridge of macroscopic size 
with pressure P = and ho = a (see ref. [16] for a separate 
derivation of this stability of the wedge considered on its 
own). The scaling also complies with the intuitive picture 
of the Rayleigh-Plateau instability as the only mechanism 
through which a macroscopic ridge breaks up into a chain 
of droplets, the size and spacing of the drops also being 
proportional to the width of the ridge. In the following we 
shall demonstrate that drops are not unstable with respect 
to deformations affecting the circular shape of the CL. 

Stability of axially symmetric drops. For the 

discussion of axisymmetric drops we introduce polar co- 



ordinates (r, ip). With a symmetric ansatz h(r) for the 
stationary solution eq. (10) turns into 

$' (h) - cr^dr [r d r h(r)] = P . (20) 

Similarly to the y-invariant case, eq. (20) admits solutions 
h(r) that have a single apex at r = (i.e., d r h(r) = 
for r = 0) and tend to a thickness ho at r = oo with 
d r h(r) < for r > 0. These are the equilibrium droplet 
shapes on which we base our stability analysis. 

In view of the rotational symmetry of the basic solution 
of eq. (20) we again focus on the Fourier components of 
the perturbations: 

*(r, ip; m) — tp(r; m) cos [m <p + a(m)] . (21) 

Equation (21) is similar to eq. (17), with the important 
difference that the continuous wave vector k is replaced by 
an integer m (resembling an azimuthal quantum number) . 

We restrict our analysis to positive m, because the de- 
formations associated with m = do not affect the circular 
shape of the contact line. For m > 1, the volume conser- 
vation constraint in eq. (9) is automatically satisfied, and 
the linear stability analysis reduces to the following ID 
eigenvalue problem for ip (r) : 

H m iP{r) = EiP{r) , H m =f + U m , (22) 

Um = U(r) + a^ T , U(r)=&'[h(r)] , (23) 

f = -<7-d r (rd r ). (24) 
r 

Instability is again indicated by the occurrence of neg- 
ative eigenvalues of H m . Before we address the overall 
stability of a stationary profile h, it is useful to make the 
following remarks about Pt m : 

1. With respect to the scalar product (ipip) = J rdrip*(p, 
the operator T is Hermitian and positive definite. 

2. For every value of m, H m is a Hermitian ID operator 
consisting of a potential term U m and a kinetic term 
T . The existence of a non-degenerate ground state 
with energy E m and without nodes (see, e.g., § 20 in 
ref. [24] and chap. 6 in ref. [25]) holds for H m . 

3. E m < E n if m < n. Indeed, for any ip (r), one has 

(H n )^ = (H m )^ +a(n 2 - m 2 ) (r~ 2 ) ^ . (25) 

Specifically, for m < n and ip chosen to be the ground 
state of H n , one has E m < (^H m ^^ < (h^j = E n . 

As in the case of the ridge, the drop is marginally stable 
with respect to lateral shifts: 

ip! (r) = d r h(r) => Hiip! = . (26) 
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The eigenfunction ipi has no nodes other than the one 
at r = 0, and thus it is the ground state of Hi; hence 
(see remark 3) E m > for all m > 1. This allows us 
to affirm that axially symmetric drops are not subject to 
instabilities which are associated with corrugations of the 
contact line. 

Conclusion. — Within the framework of a semi- 
microscopic interface displacement model we have ana- 
lyzed the linear stability of stationary distributions of a 
non-volatile liquid on a homogeneous, partially wet sub- 
strate. We have examined macroscopically extended fila- 
ments (ridges) and axially symmetric drops. The stability 
analysis has been focused on deformations affecting the 
shape of the CL which without perturbation is straight for 
ridges and circular for drops. Within our approach the CL 
tension is borne out as a microscopic property associated 
with the intrinsic morphology of the three-phase contact 
region, rather than being an arbitrary macroscopic input 
parameter. 

For each of the two aforementioned geometries, we have 
considered the elementary perturbations which induce pe- 
riodic deformations of the CL. We have shown that these 
perturbations correspond to a positive second variation 
of the free energy of the liquid structures, with the sole 
exception of the expected Rayleigh-Plateau instability for 
ridges. Notably, the critical wavelength of the Rayleigh- 
Plateau instability is not affected by CL tension. There- 
fore the interface displacement model leaves no room for 
short-wavelength CL instabilities. This stability holds for 
generic liquid-substrate and liquid-liquid interactions cap- 
tured by the effective interface potential $, which allows 
for arbitrary values and either sign of the CL tension. This 
approach also covers the limiting cases for which the CL 
tension is ill-defined according to eq. (7): e.g., for a film 
of nematic liquid crystals with antagonistic anchoring $ 
includes a long-ranged elastic term so that $(/i) <~ 1/h for 
h — > oo (see, e.g., ref. [28]). 

We conclude that the short wavelength instability found 
in ref. [20] for wavelengths shorter than the microscopic 
width of the CL region is outside the range of the macro- 
scopic capillary model employed, and thus not physical. 
Our analysis within the framework of the interface dis- 
placement model generalizes the result of ref. [16] and 
shows that, at least for small contact angles, negative CL 
tensions do not lead to short wavelength instabilities. 

* * * 
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